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Focusing on a quantum-limit behavior, we study a sin- 
gle vortex in a clean s-wave type-II superconductor by self- 
consistently solving the Bogoliubov-de Gennes equation. The 
discrete energy levels of the vortex bound states in the quan- 
tum limit is discussed. The vortex core radius shrinks mono- 
tonically up to an atomic-scale length on lowering the tem- 
perature T, and the shrinkage stops to saturate at a lower T. 
The pair potential, supercurrent, and local density of states 
around the vortex exhibit Friedel-like oscillations. The local 
density of states has particle-hole asymmetry induced by the 
vortex. These are potentially observed directly by STM. 

PACS number(s): 74.60. Ec, 61.16.Ch, 74.72.-h 

Growing interest has been focused on vortices both in 
conventional and unconventional superconductors from 
fundamental and applied physics points of view. This is 
particularly true for high-T c cuprates, since it is essen- 
tial that one understands fundamental physical proper- 
ties of the vortices in the compounds to better control 
various superconducting characteristics of some techno- 
logical importance. Owing to the experimental devel- 
opments, it is not difficult to reach low temperatures 
of interest where distinctive quantum effects associated 
with the discretized energy levels of the vortex bound 
states are expected to emerge. The quantum limit is 
realized at the temperature where the thermal smear- 
ing is narrower than the discrete bound state levels |Q: 
T/T c < 1/(/cf£o) with £o=vf/A the coherence length 
(Ao the gap at T — 0) and k-p (vf) the Fermi wave num- 
ber (velocity). For example, in a typical layered type-II 
superconductor NbSe2 with T c = 7.2 K and /cf£o ~ 70, 
the quantum limit is reached below T < 50 mK. As for 
the high-X^ cuprates, the corresponding temperature is 
rather high: T < 10 K for YBa 2 Cu307-,5 (YBCO). 

Important microscopic works to theoretically investi- 
gate the quasiparticle spectral structure around a vortex 
in a clean limit are put forth by Caroli et al. ^ , Kramer 
and Pesch Q|, and Gygi and Schliiter ||. The low-lying 
excitations are essential to correctly describe low-T ther- 
modynamic and transport properties in the vortex state 
(or the mixed state). These include anomalous electric 
or thermal Hall conductivity Q and mysterious obser- 
vations of the quantum magnetic dHvA oscillations [^J; 
various topics are debated intensively M. Yet there has 
been no serious attempt or quantitative calculation to 
explore deep into the quantum regime. 

The purposes of the present paper are to reveal the 
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quantum-limit aspects of the single vortex in s-wave su- 
perconductors and to discuss a possibility for the obser- 
vation of them. 

The present study is motivated by the following re- 
cent experimental and theoretical situations: (1) The so- 
called Kramer and Pesch (KP) effect [0J|J^,|]; a shrink- 
age of the core radius upon lowering T (to be exact, an 
anomalous increase in the slope of the pair potential at 
the vortex center at low T) is now supported by some ex- 
periments U . The T dependence of the core size is stud- 
ied by /j,SR on NbSe2 and YBCO jl(|, which is discussed 
later. The KP effect, if confirmed, forces us radically al- 
ter the traditional picture jllj for the vortex line such 
as a rigid normal cylindrical rod with the radius £o- (2) 
The scanning tunneling microscopy (STM) experiment 
on YBCO by Maggio-Aprile et al. |L2|], which enables us 
to directly see the spatial structure of the low- lying quasi- 
particle excitations around the vortex, arouses much in- 
terest. They claim that surprisingly enough, there exist 
only a few discretized bound-state levels in the vortex 
core, i.e., the vortex is almost "empty." It resembles 
our naive image for conventional s-wave superconduc- 
tors where in the quantum limit a few quantized levels 
of the bound states remain inside the bulk energy gap 
Ao- (3) The theoretical situation on this subject Jl2j is 
still very confusing; Some Jl3| claim that the bound-state 
energy levels are not discretized for e?-wave pair, but dis- 
cretized for s-wave pair. Some |l4| claim the discretized- 
like structure even for the former. For s-wave case, where 
the formulation of the problem is well defined, we should 
establish our understanding of the vortex structure in the 
quantum limit. (4) Lastly, we are motivated by a curios- 
ity; Previously we have calculated the local density of 
states (LDOS) for s-wave pair on the basis of the qua- 
siclassical (Eilenberger) theory |^5| , successfully applied 
to the STM observations on NbSe2 done by Hess et al. 
p6fl . We are particularly interested in what happens in 
LDOS at further lower T, say, below 50 mK deep into 
the quantum limit, at which it may be now feasible to 
perform STM experiments. 

Prompted by these motivations, we self-consistently 
solve the Bogoliubov-de Gennes (BdG) equation, which 
is one of the most fundamental microscopic equations 
of superconductivity and contains fully quantum effects. 
We start with the BdG equations for the quasiparticle 
wave functions Uj(r) and Vj(r) labeled by the quantum 
number j: 
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Uj (r) + A(r)uj(r) = E jUj (r), 
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v j (T) + A*(r)u j {r)=E j v j (r), (1) 



in a dimensionless form, where A(r) is the pair potential 
and Ep (—kp^ /2) the Fermi energy. The length (energy) 
scale is measured by £ (Ao)- For an isolated single vortex 
in an extreme type-II superconductor, we may neglect 
the vector potential in Eq. ([!]). The pair potential is 
determined self-consistently by 



A(r)=g «i(r)t$(r){l- 2/(25,-)} 



(2) 
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with the Fermi function f(E). Here, g is the coupling 
constant and wd the energy cutoff, which are related by 
a BCS relation via the transition temperature T c and 
the gap A . We set cup, = 10A . The current den- 
sity is given by j(r) oc lmJ2j [/(i?j)uJ(r)V%(r) + {1 - 
f(Ej)}vj(r)\7v*(r)]. We consider an isolated vortex un- 
der the following conditions, (a) The system is a cylinder 
with a radius R. (b) The Fermi surface is cylindrical, 
appropriate for the materials such as NbSe2 and high-T c 
cuprates. (c) The pairing has isotropic s-wave symmetry. 
Thus the system has a cylindrical symmetry. We write 
the eigenf unctions as Uj(r) = u nfi (r) exp\i(fj> — \)0\ and 
Vj(r) = u„ M (r)exp[i(/i+!)6>] with A(r) = A(r) exp[-£0] 
in polar coordinates, where n is a radial quantum num- 
ber and the angular momentum = 3, 1, We 
expand the eigenfunctions in terms of the Bessel func- 
tions J m (r) as u nil (r) = J2i c m<t>i\v,-±\{ r ) and = 
E 4 d ™^|M+i|( r ) with 4> lm {r) = RJin ^ ia . m) Jm{a im r/R), 
(i = 1, 2, • • • , N, and cn m is the i-th zero of J m (r)) ■ The 
BdG is reduced to a 2N x 2N matrix eigenvalue problem 
pL Our system is characterized by /cf£oj which is a key 
parameter of the present problem. 

In Fig. [I], the calculated spatial variation of A(r) is 
displayed for various T. It is seen that as T decreases, 
the core size £1 defined by ^ x = lim r ^ A (r)/r shrinks 
and the oscillatory spatial variation with a wave length 
~ 1/fcp becomes evident in A(r) The physical rea- 

son for this Friedel-like oscillation lies in the following 
facts. All eigenfunctions u n fj,(r) and v niJ ,{r) contain a 
rapid oscillation component with 1/fcp. At lower T the 
lowest bound states, whose oscillation amplitude is large 
near the core, dominate physical quantities. We note that 
the oscillatory behavior can always appears at sufficiently 
low T irrespective of values of kp^Q. We also mention that 
a similar oscillatory spatial variation around a vortex core 
in the Bose condensate of 4 He is found theoretically, due 
to the roton excitations jlTj . 

The associated supercurrent jg(r) and the field H(r) 
are shown in Fig. 0. Reflecting the above oscillation, jo (r) 



also exhibits a weak oscillation around r — 0.2-0.5 £o- It 
is difficult to see the oscillation in H(r), because it is 
obtained by integrating j$(r) via the Maxwell equation 
V x H = ^j(r), resulting in a smeared profile. It is 
also seen that the position of the maximum of jg (r) be- 
comes shorter as T decreases. These features quite differ 
from those obtained within the Ginzburg-Landau (GL) 
framework 

The T dependence of £i(T) for various kp^o values is 
shown in Fig. |3[ Coinciding with Kramer and Pesch 
|l| for s-wave pair and Ichioka et al. || for d-wave 
pair, £i(T) decreases almost linearly with T, that is, 
£i( r )/£o ~ T/T c except at extremely low T. An im- 
portant difference from these quasiclassical theories 
appears at lower T. At a lower T < T ~ T c /(kp£ ), 
where the quantum limit is realized, the shrinkage of the 
core size stops to saturate, and the saturated value is 
estimated as £i/£o ~ (^fCo) -1 - 

According to the /iSR experimental data jl(j , the core 
radius in NbSe2 shows a strong T dependence, while that 
in YBCO with T c =60 K is almost T-independent below 
^0.6T C . This seemingly contradicting result can be un- 
derstood as follows. The strong T dependence in NbSe2 
is the usual KP effect corresponding to the curves for 
larger /cf£o in Fig. [| At lower T than To estimated 
as ^100 mK (fcp£o ~ 70), the shrinkage must saturate 
(the above experiment is done above ~2 K). As for the 
YBCO data, since the estimated kp^Q is small (~4 E^] 
for YBCO with T c =90 K), the saturation is already at- 
tained at a relatively high T such as shown in Fig. |^. 
Thus the absence or weakness of the KP effect in YBCO 
is simply attributable to the fact that the quantum-limit 
temperature To is quite high. 

Reflecting the shrinkage of the core radius, the bound- 
state energies E^ increases as T decreases. This In- 
dependent E^ shift, due to the KP effect, and its satu- 
ration at lower T may lead to a nontrivial T dependence 
in thermodynamic and transport properties. 

In Fig. [|, we plot the energy levels E^ of the low- 
lying bound states hu = |> §)"*'> if) as a function of 
fcp£oi a t sufficiently low T (T/T c = 0.01) where increas- 
ing of the energy levels saturates. It is seen that in 
large-fcp£o region, the bound states densely pack inside 
the gap Aq, allowing us to regard them as continuous 
ones. This is the case where the quasiclassical approx- 
imation validates. In small-fcp^o region, where the 
quantum effect is important even at high T, only a few 
bound states remain within the low-energy region. We 
find that even in small- \fi\ region, the spacing between 
the energy levels E^ is not constant, but rather becomes 
narrower as increases. The often adopted formula 
E^/Aq = 2^/(fc F £o) or 2/z/(fc F £i) due to Caroli et al. M, 
or E^/Aq — (2(i/kp£o) ln[£o/2£i] by Kramer and Pesch 
in the limit £1 <c £0 jy do not satisfactorily explain our 
self-consistent results. Instead, our result is empirically 
fitted to a formula Ei/ 2 /Aq = (0.5/&f£o) ln[fcp£o/0.3] for 
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large fcpCo as shown in the dotted curve in Fig. |[ 

In Fig. [|, the spectral evolution, i.e., the spatial 
variation of LDOS, which is calculated by N(r, E) oc 
E,-[k-(r)| 2 /'(£ - E S ) + Mr)| 2 /'(£ + £,•)], * shown 
for fcp£o=8 at low temperature T=0.05T C . It is well 
contrasted with that of the higher T case by Gygi and 
Schliiter || (see, for comparison, Fig. 15 in Ref. ||] 
where /efCo ~ 70 and T ~ 0.13T C , calculated under 
the two-dimensional Fermi surface). As lowering T, be- 
cause of the quantum effects, the thermally smeared spec- 
tral structure drastically changes and becomes far finer 
one around the vortex. The spectra are discretized in- 
side the gap and consist of several isolated peaks, each 
of which precisely corresponds to the bound states 
(|/i| — ^, |, • • •)• Reflecting the oscillatory nature of the 
eigenf unctions u^ir) and v^if) with the period the 
spectral evolution also exhibits the Friedel-like oscillation 
as seen from Fig. |[ 

To show clearly the particle-hole asymmetry of the 
LDOS of Fig. J| which is another salient feature, we 
present in Fig. g the spectra at the vortex center r = 
and 0.2£o [We can barely see the asymmetry in Wang and 
MacDonald |@ (see Fig. 3(a) in Ref. Q)] . At the cen- 
ter r = 0, the bound-state peak with £1/2 appears only 
at E > side, because the eigenfunction U\/i{r = 0) 7^ 0, 
which consists of the Bessel function Jo(r = 0) (^ 0), and 
all others for \E^\ < Ao vanish at r — 0. At r = 0.2Co, 
the other bound-state peaks are seen. The particle-hole 
asymmetry in the vortex bound states appears even if 
the normal-state density of states is symmetric. These 
features are subtle || or absent |15| in the previous cal- 
culations. This asymmetry around the vortex is quite 
distinctive, should be checked by STM experiments, and 
may be crucial for the Hall conductivity in the mixed 
state. 

Let us argue some of the available experimental data 
in the light of the present study. The lowest bound state 
level E1/2/A.0 is estimated by Maggio-Aprile et al. |ll| for 
YBCO with T c =90 K (E 1/2 =5.5 meV and A =20 mcV), 
yielding /cfCo ~ 4. Since it implies that £0 is only of 
the order of the crystal-lattice constant, we should cau- 
tion that Maggio-Aprile et al. p2] take their data for 
the spectral evolution every 10 A apart near the core, 
thus the important spatial information on LDOS might 
be lost. So far the existing STM data fl|,|l6],[l9) taken 
at the vortex center are almost symmetric about E=0, 
e.g., on NbSe2 at T=50 mK |jl6). The reason why the 
so-called zero-bias peak is centered just symmetrically at 
E=Q is that A;fCo is large and T is too high to observe 
the quantum effects. 

We emphasize that in any clean s-wave type-II super- 
conductors at appropriately low T (< To ~ T c /(fcpCo)) 5 
one can observe these eminent characteristics associated 
with the quantum effects. For example, a typical A-15 
compound V 3 Si [§ with T c =17 K, Co ^ 60 A (fc F Co ^ 



12), and a borocarbide LuNi 2 B 2 C [|J| with T c =16 K, 
Co — 80 A are the best candidates to check our results. 

In summary, we have analyzed the vortex core struc- 
ture and the related quasiparticle energy spectrum by 
self-consistently solving the BdG equation for an isolated 
vortex in a clean s-wave type-II superconductor, focus- 
ing on the low-T quantum effects. We have found the far 
richer structure in the pair potential, supercurrent, and 
LDOS than what one naively imagines from the corre- 
sponding calculations done at high T or fc F Co 3> 1 ||,[l5) > 
and pointed out experimental feasibility to observe it. 

The widely used working hypothesis for the vortex core 
of a rigid normal rod with the radius Co jll| must be 
cautiously used for the clean superconductors of inter- 
est: the magnetic field distribution probed by neutron 
diffraction ^FJ or /iSR jl0| through the magnetic form 
factor analysis based on the GL theory must be taken 
with caution. Detailed investigations of various myster- 
ies associated with the vortices, e.g., the thermal Hall 
conductance PJ belong to future work. 

We would like to thank J. E. Sonier and A. Yaouanc 
for useful discussions. 

FIG. 1. The spatial variation of the pair potential A(r) 
normalized by Ao around the vortex for several temperatures 
and /cfCo = 16. The length r is measured by £o- 



FIG. 2. The current distribution normalized by 

c^>o/(87r 2 £oK 2 ) for several temperatures, where cj>o is the flux 
quantum and k (3>1) is the GL parameter. The inset shows 
the field distribution normalized by 0o/(27t£qK 2 ). The tem- 
peratures are the same as in Fig. 1, and /cf^o = 16. 



FIG. 3. The T dependence of the vortex radius £1 nor- 
malized by £0 for several ^fCo (= 1.2, 2, 4, and 16 from top 
to bottom). 



FIG. 4. The lowest seven bound-state energies £7 M , nor- 
malized by Ao, as a function of kp^o, at enough low temper- 
ature T/T c = 0.01. The dotted line is a fitting curve (see the 
text). 



FIG. 5. The spectral evolution N(E,r) at T/T c = 0.05 
and fcF^o = 8. It is normalized by the normal-state density 
of states at the Fermi surface. E and r are measured by Ao 
and Co, respectively. 



FIG. 6. The local density of states N(E, r) at r = (solid 
line) and 0.2Co (dotted line). T/T c = 0.05 and fc F Co = 8. 
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